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Abstract 

Hypersonic  endoatmospheric  interceptors 
are  presently  under  serious  consideration  for  stra¬ 
tegic  and  theater  missile  defenses.  On  board  infra¬ 
red  seekers  are  planned  for  terminal  guidance.  It  is 
desirable  to  perform  seeker  trade  studies  determin¬ 
ing  such  things  as  expected  S/N,  tracking  perfor¬ 
mance  and  optimum  wavelength  of  operation.  The 
dynamics  of  the  interceptor/ target  trajectories  com¬ 
bine  with  atmospheric  interaction  at  both  vehicles  to 
degrade  the  infrared  signal  of  the  target.  Because  of 
this  degradation,  computation  of  this  signal  at  the 
seeker  focal  plane  and,  consequently,  performance  of 
the  seeker  trade  studiesare  made  very  difficult.  These 
computations  are  usually  performed  piecemeal,  fora 
limited  number  of  trajectory  points,  and  for  very 
specific  seeker  configurations.  An  integrated  radi¬ 
ance  throughput  model  is  under  development  that 
encompasses  the  trajectory  dynamics  and  the  atmo¬ 
spheric  interactions.  The  model  is  intended  for 
medium  fidelity  performance  analysis  to  enable  over¬ 
all  systems  performance  evaluations  and  trades.  Per¬ 
formance  of  proposed  window  materials  under  a 
variety  of  aerothermal  environments  is  an  example  of 
a  typical  trade  space.  Various  existing  codes  are 
utilized  integrally  within  the  model,  or  are  used  to 
generate  lookup  tables  for  parameters  over  the  range 
of  trajectory  and  atmospheric  variables.  These  in¬ 
clude:  the  Optical  Signature  Code  and  the  SIRRM 
Code  for  target  signatures  and  bow  shock  radiance. 
The  LOWTRAN  7  Code  is 

used  integrally  to  calculate  atmospheric  transmis¬ 
sion  and  background.  S/N  and  tracking  perfor¬ 
mance  will  be  estimated  using  a  simple  detector 
module.  The  model  is  described,  and  initial  results 
are  given  that  illustrate  atmospheric  transmission 
and  window  performance  effects. 


Introduction 

A  re-entering  RV  rapidly  heats  up  compris¬ 
ing  a  dynamic  target  signature.  The  target  signature 
propagates  through  a  changing  atmospheric  path 
(with  subsequent  loss  of  energy  through  scattering 
and  absorption)  arriving  at  the  hypersonic  intercep¬ 
tor.  This  hypersonic  interceptor  also  encounters  a 
dynamic  and  extreme  aerothermal  environment.  The 
interceptor  seeker  window  heats  up,  and  may  re¬ 
quire  cooling  by  flowing  cold  gas  over  its  surface. 
Simultaneously,  the  bow  shock  and  the  heated  win¬ 
dow  of  the  interceptor  radiate  and  produce  back¬ 
ground  noise.  All  of  these  effects  contribute  to  deg¬ 
radation  of  the  target  image  in  the  seeker  focal  plane. 

The  dynamics  described  above  comprise  only 
a  fraction  of  the  effects  that  must  be  considered  in 
evaluating  the  performance  of  a  seeker  on  board  an 
endoatmospheric  hypersonic  interceptor.  Tradition¬ 
ally,  the  problem  of  determining  the  signal  and  back¬ 
ground  incident  on  a  detector  is  done  in  piecemeal 
fashion,  and  often  does  not  encompass  all  of  the 
signal  degrading  factors.  Furthermore,  new  technol¬ 
ogy  development,  ongoing  research  and  test  pro¬ 
grams  are  shedding  new  light  on  many  of  the  effects 
encountered  when  imaging  from  a  hypersonic  ve¬ 
hicle.  Therefore,  a  need  exists  for  a  model  which  will 
examine  all  of  the  effects  that  degrade  the  interceptor 
seeker  performance,  and  that  can  do  so  for  a  wide 
range  of  environments,  and  for  a  variety  of  seeker 
technologies  and  concepts. 

The  Integrated  Radiance  Throughput  Model 
(IRTM)  is  being  developed  to  quickly  and  efficiently 
provide  an  integrated  radiance  throughput  analysis, 
considering  all  sources  of  signal  degradation,  for 
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Figure  1  -  Target  To  Detector  Transmission 


endoatmospheric  interceptor  vehicle  seekers.  The 
model  will  be  used  to  examine  the  wavelength  of 
operation  of  electro-optical  (EO)  seekers  as  new  win¬ 
dow  and  seeker  technologies  are  being  developed. 
Simultaneously,  the  model  will  incorporate  data  gen¬ 
erated  by  new  and  innovative  CFD  and  wave  optic 
codes.  Data  from  tests  (such  as  the  upcoming  aero- 
optics  tests  at  Calspan)  can  be  used  as  inputs,  and  can 
help  validate  the  model.  The  model  is  modular  in 
nature,  and  includes  options  within  various  modules 
to  provide  analyses  over  a  wide  range  of  conditions, 
inputs  and  levels  of  fidelity. 

The  initial  thrust  of  the  program  was  to  pro¬ 
vide  a  means  of  evaluating  diamond  window  and 
seeker  technologies  being  developed  under  a  Strate¬ 
gic  Defense  Command  endo- technology  effort.  It  has 
subsequently  been  expanded  to  include  a  wide  range 
of  window  materials  and  technologies.  The  earliest 
version  provides  a  basic  throughput  of  the  signal 
from  a  target,  through  the  atmosphere,  and  through 
a  window.  It  is  presently  being  used  to  compare 
performance  of  various  window  materials  over  a 
range  of  target  and  atmospheric  path  parameters. 

Later  versions  will  include  additional  mod¬ 


ules  (i.e.  consider  more  signal  degrading  sources), 
will  have  higher  fidelity,  and  will  calculate  through¬ 
put  for  multiple  points  within  the  relative  intercep¬ 
tor/target  trajectories.  The  final  version  will  have 
provisions  for  inputting  seeker  and  aero-optic  pa¬ 
rameters  derived  from  analysis  or  testing  of  actual 
hardware.  The  intent  of  this  paper  is  to  describe  the 
overall  model,  discuss  the  present  statusof  the  model, 
outline  the  plans  for  completing  the  various  mod¬ 
ules,  and  to  display  initial  results  in  which  several 
window  materials  are  compared. 


Integrated  Radiance  Throughput  Model  Description 

Figure  1  indicates  a  general  picture  of  the 
target  signature  and  sources  of  image  degradation 
that  are  being  incorporated  in  the  model.  As  the 
target  progresses  through  its  trajectory  it  rapidly 
heats,  and  displays  a  distinct  radiative  signature  that 
is  primarily  dependent  on  the  temperature.  The 
model  can  include  background  created  from  nuclear 
events  as  well  as  bow  shock  radiance  background. 
Both  of  these  are  calculated  by  the  SIRRM  Code  and 
areinputinto  the  model  as  lookup  tables.  LOWTRAN 
is  integral  to  the  model  and  is  used  to  calculate 
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Figure  2  -  Model  Module  Breakdown 


atmospheric  effects.  Initially,  very  simple  models 
will  be  used  to  simulate  the  aero-optics  effects  of  the 
turbulent  coolant  flow.  The  inset  to  figure  1  illus¬ 
trates  multiple  image  degrading  effects  that  take 
place  at  the  interceptor. 

Figure  2  illustrates  all  of  the  modules  that 
will  be  included  in  the  completed  model.  Note  the 
several  modules  that  are  required  to  account  for 
effects  right  at  the  interceptor.  The  window  module 
must  account  for  change  in  absorption  as  a  function 
of  heating,  background  caused  by  self-emission  of 
the  window,  and  distortion  in  the  window.  Bow 
shock  radiance  and  aero-optics  effects  are  accounted 
forby  a  combination  of  three  modules:  1 )  the  aerody¬ 
namic  and  coolant  flow  effects  module,  2)  the  bow 
shock  radiance  module,  and  3)  the  aero-optic  image 
degradation  and  jitter  module. 

The  modules  will  consist  of  either  1 )  industry 
standard  computer  codes,  2)  data  bases  generated  by 
computer  codes  or  actual  test  data,  or  3)  algorithms 
based  on  the  physics  of  the  particular  module.  As  an 
example:  LOWTRAN  is  used  to  directly  calculate 
atmospheric  transmission.  However,  it  is  also  used 
to  genera  tedata  tables  for  representative  atmospheric 
and  earthshine  backgrounds.  It  is  necessary  to  do  the 


latter  as  lookup  tables  because  of  very  long  run  times 
of  the  LOWTRAN  program  when  calculating  back¬ 
ground.  It  is  desired  that  the  model  operate  rela¬ 
tively  quickly  so  that  multiple  points  in  a  trajectory 
can  be  evaluated,  and  multiple  technologies  can  be 
evaluated  comparatively. 

The  target  signal  mod  ule  is  another  example. 
Initial  cases  are  run  using  Planck's  black  body  equa¬ 
tions  to  calculate  radiance  at  a  given  target  tempera¬ 
ture.  Later,  higher  fidelity  models  will  use  the  Opti¬ 
cal  Signature  Code  to  generate  look  up  tables  that 
more  accurately  represent  target  signatures  under 
various  velocity  and  altitude  conditions. 

The  mission  module  serves  as  an  input  and 
coordination  module.  The  multipoint  version  of  the 
model  will  have  relative  trajectories  input  into  the 
mission  module.  In  turn,  it  will  calculate  altitudes, 
velocities  and  relative  ranges  between  the  intercep¬ 
tor  and  target.  The  RV  altitude  and  velocity  param¬ 
eters  will  be  used  to  access  appropriate  target  signa¬ 
tures.  The  relative  range  parameters  will  feed  the 
LOWTRAN  program  to  enable  appropriate  trans¬ 
mission  calculations. 
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Figure  3  -  Sample  Throughput 


Aerodynamic  and  coolant  flow  effects  will 
be  calculated  off-line  by  CFD  codes  to  establish  in¬ 
puts  into  the  SIRRM  code.  The  SIRRM  code  will  also 
be  run  off-line  to  generate  lookup  tables  that  are 
accessed  by  using  appropriate  interceptor  altitude 
and  velocity  parameters.  CFD  code  outputs  may  also 
be  used  to  generate  inputs  for  various  aero-optic 
codes.  Aero-optic  degradations  indicated  by  these 
codes  can  be  input  into  the  throughput  model  as 
simple  PSF  degradation  factors. 

The  detector  module  will  consist  of  either  a 
quadrant  detector  model  with  appropriate  spectral 
andresponsivityparameters,ormaybeimplemented 
as  a  small  staring  array  model  with  a  simple 
centroiding  algorithm. 


Window  distortion  can  be  input  as  look  up 
tables  generated  by  using  thermal  modelling  codes 
that  in  turn  are  fed  appropriate  temperature  input 
parameters  by  CFD  code  results. 

In  general,  industry  standard  codes  that  run 
sufficiently  fast  will  be  utilized  directly  in  the  model. 
Large  codes  that  require  extensive  computational 


power  will  be  used  offline  to  generate  lookup  tables. 
This  way,  the  overall  model  can  operate  rapidly  and 
evaluate  performance  over  a  range  of  mission  param¬ 
eters,  and  can  be  used  to  evaluate  several  seeker 
technology  combinations.  The  model  is  designed  to 
operateonPCs.  A  Macintosh  version  is  also  planned. 

Present  Status  and  Module  Implementation 

The  present  version  of  the  model  is  imple¬ 
mented  to  calculate  signal  throughput  from  a  target, 
through  the  atmosphere,  and  through  a  seeker  win¬ 
dow.  The  intent  of  this  version  is  to  compare  perfor¬ 
mance  of  various  window  materials  given  specific 
target  signatures  and  atmospheric  path  geometries 
and  conditions,  and  assumed  window  temperatures. 

Figure  3  indicates  the  process  involved  in 
this  simplified  analysis.  The  target  is  modeled  as  a 
black  body.  The  user  chooses  the  desired  target 
temperature,  the  RV  and  interceptor  altitudes  and 
the  slant  range.  The  desired  atmospheric  model  is 
chosen  within  the  LOWTRAN  suite  of  options.  The 
altitude  and  slant  range  information  is  used  with 
LOWTRAN  to  calculate  an  atmospheric  transmit¬ 
tance.  The  target  signature  undergoes  a  spectral 
multiplication  by  the  calculated  atmospheric  trans- 
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mittance  factor.  It  is  further  multiplied  by  the  1  /R 
factor  to  determine  signal  irradiance  at  the  intercep¬ 
tor.  Transmittance  through  the  window  is  calculated 
and  multiplies  this  signature.  The  resulting  signa¬ 
ture  is  then  multiplied  by  the  aperture  area  and  the 
transmission  factor  of  the  optics.  The  result  is  a 
calculation  of  signal  intensities  resident  on  the  detec¬ 
tor. 

Plots  of  the  result  from  each  step  can  be 
produced  using  the  data  tables  generated  by  the 
model.  Figure  3  specifically  depicts  a  calculation  for 
a  hypothetical  diamond  window. 

Figure  4  depicts  the  flow  diagram  for  those 
modules  implemented .  As  can  be  seen  this  is  a  subset 
of  the  full  model  illustrated  in  Figure  2.  A  final  note 
on  present  status:  initial  window  property  calcula¬ 
tions  used  lookup  tables  of  absorption  properties  of 
materials  found  in  literature  or  in  various  govern¬ 
ment  contract  reviews.  We  recently  obtained  a  copy 
of  the  "PhononB"  program  [Mike  Thomas,  Johns 
Hopkins  Applied  Physics  Laboratory],  and  will  use 
it  to  generate  new  lookup  tables. 

The  target  module  currently  calculates  sig¬ 
nature  using  Planck's  Black  Body  Law, 

WL=  (K1  x  AXx  AT)/(Xs(eK2AT-l))  0) 

where  K1  =  37405,  K2  =  14378.9,  A,  is  the  radiating 
area  incm2,X  is  wavelength  in  pm,  T  is  temperature 
in  K°  and  WL  is  power  radiated.  The  power  radiated 
is  calculated  at  each  wavelength  increment  over  the 


total  waveband  (which  is  specified  as  an  input). 

The  atmospheric  attenuation  and  back¬ 
ground  module  uses  the  LOWTRAN  computer  code 
to  determine  the  atmospheric  transmission  of  the 
target  signal  over  the  same  waveband  and  wave¬ 
length  increment  as  specified  in  the  target  module. 
The  interceptor  altitude,  target  altitude  and  slant 
range  are  input  using  the  LOWTRAN  input  screen. 
The  atmospheric  transmission  values  calculated  by 
LOWTRAN  are  then  multiplied  by  the  power  radi¬ 
ated,  thereby  calculating  the  amount  of  signal  trans¬ 
mitted  through  the  atmosphere.  This  value  is  also 
multiplied  by  the  1  /RJ  factor  to  account  for  the  range 
between  the  target  and  the  interceptor  (eqn.  2). 

Wt=  (WLx  ta)/R2  (2) 

where  W(  is  the  power  incident  at  the  interceptor,  ta  is 
the  atmospheric  transmission  value,  and  R  is  the 
range  between  the  target  and  the  interceptor. 

The  window  module  contains  absorption 
data  (vs.  wavelength)  for  several  different  materials 
including  diamond  and  sapphire,  two  materials  be¬ 
ing  used  for  endoatmospheric  interceptors.  After 
choosing  a  window  material  and  window  thickness, 
the  power  transmitted  through  the  window  is  calcu¬ 
lated  by 

Ww  =  W1xe-(*X,)  <3) 

where  Ww  is  the  power  transmitted  through  the 
window,  a  is  the  absorption  coefficient  of  the  win- 


Power 


Figure  4  -  Throughput  Model  Flow 
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dow  material,  and  t  is  the  thickness  of  the  window. 

The  final  module  implemented  in  the 
throughput  model,  is  the  detector /optics  module.  It 
currently  calculates  the  power  received  per  quadrant 
of  the  detector  by 


trade  factors  may  be  more  important,  i.e.,  much  su¬ 
perior  material  properties  for  diamond  that  may 
simplify  other  aspects  of  the  overall  interceptor/ 
seeker  design  (diamond  is  stronger,  much  better 
thermal  conductance,  has  lower  emissivity,  etc.). 


WD  =  ArxtD  x  Ww/4  (4) 

where  Ar  is  the  area  of  the  receiver  and  ^  is  the 
transmission  factor  of  the  optics.  For  the  results 
shown  in  the  following  paragraph,  the  aperature 
diameter  is  10cm  and  the  transmission  of  the  optics 
is  assumed  to  be  .8. 

Figures  5-7  show  a  comparison  of  diamond 
and  sapphire  windows  looking  at  a  900°K  target 
from  50km  away.  The  sapphire  window  represents 
a  typical  HEDI  configuration  and  the  diamond  win¬ 
dow  represents  a  typical  Endo  LEAP  configuration. 
Figure  5  shows  the  results  for  a  target  at  37km  alti¬ 
tude  and  an  interceptor  at  7km  altitude.  As  is  ex¬ 
pected,  the  sapphire  window  outperforms  the  dia¬ 
mond  in  the  waveband  from  about  3pm  -  5pm,  but 
diamond  performs  much  better  across  the  full  spec¬ 
trum,  excluding  4pm  -  6pm.  Figure  6  shows  the 
results  for  the  target  at  40km  altitude  and  the  inter¬ 
ceptor  at  10km  altitude.  The  irradiance  for  both 
window  ma  terials  has  improved  considerally  a  t  these 
higher  altitudes  (because  of  the  less  dense  atmo¬ 
sphere  to  transmit  through),  most  notably  for  the 
diamond  at  the  longer  wavelengths.  Figure  7  shows 
the  diamond  and  sapphire  windows  for  the  target  at 
45km  and  the  interceptor  at  15km.  Once  again,  the 
diamond  performance  is  greatly  increased  at  longer 
wavelengths.  Figures  8-10  show  the  power  received 
per  quadrant  of  the  detector  for  the  same  cases  out¬ 
lined  in  Figures  5-7.  Three  spectral  bands  were 
input:  band  1  from  3-5pm,  band  2  from  5-8pm  and 
band  3  from  8-12pm.  In  each  case,  the  diamond 
window  outperforms  the  sapphire  window  at  the 
higher  wavelengths,  (spectral  bands  2  and  3),  and 
the  sapphire  significantly  outperforms  the  diamond 
in  spectral  band  1. 

If  one  shifts  to  a  2-4  pm  band  as  shown  in 
Figure  11,  more  energy  is  available  using  a  diamond 
window.  However,  sapphire  still  outperforms  dia¬ 
mond,  but  the  results  are  closer.  Two  things  to  keep 
in  mind  in  making  trades  of  this  nature:  1)  for  targets 
of  these  temperatures  there  is  plenty  of  power  in 
these  bands  (2-4  &  3-5)  using  either  window,  2)  other 


Figure  5  -  7km  Interceptor  Throughput 


Figure  6  -  10km  Interceptor  Throughput 
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900  DEG.  TARGET  9  45  KM 
INTERCEPTOR  9  15  KM  (SO  KM  SLANT  RANGE ) 


Figure  7  -  15km  Interceptor  Throughput 
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Figure  8  -  7km  Interceptor  Detector  Power 
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Figure  10  -  15km  Interceptor  Detector  Power 
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2 . 4  MICRON  BAND 

Figure  11  -  Detector  Power  for  2-4  micron  Band 
Summary 

An  integrated  radiance  throughput  model 
is  under  development.  This  model  is  intended  for 
use  as  a  performance  evaulation  tool  for  various 
seeker  component  and  systems  technologies.  It  can 
be  used  as  a  systems  tool  to  provide  a  means  of 
performing  trade  studies  at  a  medium  level  of  fidel¬ 
ity.  It  is  already  being  used  for  comparative  studies 
of  possible  window  materials.  The  model  is  imple¬ 
mented  to  provide  rapid  results  for  numerous  mis¬ 
sion  scenarios,  and  geometries  using  various  emerg¬ 
ing  window  and  seeker  technologies. 
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Figure  9  -  10km  Interceptor  Detector  Power 
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